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Abstract 

The partial-wave expansion used to treat the distortion of scattered elec- 
trons by the nuclear Coulomb field is simpler and considerably less time- 
consuming when applied to the production of muons and electrons by low and 
intermediate-energy neutrinos. For angle-integrated cross sections, however, 
a modification of the "effective-momentum-transfer" approximation seems to 
work so well that for muons the full distorted-wave treatment is usually unnec- 
essary, even at kinetic energies as low as an MeV and in nuclei as heavy as lead. 
The method does not work as well for electron production at low energies, but 
there a Fermi function usually often proves perfectly adequate. Scattering of 
electron-neutrinos from muon decay on iodine and of atmospheric neutrinos 
on iron are discussed in light of these results. 
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1 Introduction 



Converting neutrinos with energies below a few hundred MeV into electrons or muons 
by scattering them from nuclei is useful in the search for new physics. Radiochemical 
detectors HI] rely on neutrino- nucleus scattering to measure the apparently too low solar 
neutrino flux. Proton-decay experiments, doubling as neutrino observatories, use neu- 
trino scattering from oxygen or iron to look for oscillations of neutrinos produced in the 
atmosphere |2], g]. And experiments at Los Alamos [Q] measure neutrino-carbon cross sec- 
tions, in part to check the flux in a beam that has reportedly produced neutrino oscillations 
in the lab. 

Because neutrinos interact so weakly the Born approximation to neutrino-nucleus scat- 
tering should be accurate. The only difficulty with this first-order plane-wave approach 
is that the electrons or muons produced when the neutrinos interact feel an electrostatic 
force from the nucleus as they leave. In heavy nuclei the Coulomb force can have large 
effects and even in light nuclei it changes cross sections noticeably. The best way to include 
the Coulomb interaction is through the distorted- wave Born approximation (DWBA)j5[|, 
which entails an expansion of the outgoing wave function in partial waves. In processes 
similar to neutrino scattering but more often measured, two approximations help avoid 
the cumbersome and time-consuming partial-wave expansion. Calculations of beta de- 
cay, which produces low-energy electrons, contain a Fermi function that multiplies the 
decay amplitude by the ratio of Coulomb to free Dirac s-wave functions at the nuclear 
radius |J. In electron scattering, at energies much larger than the electrostatic potential 
energy in the nucleus, an "effective momentum approximation" (EMA) uses plane waves 
with shortened wavelengths and increased amplitudes in place of incoming and outgoing 
distorted waves 0. Both approximations are substantial improvements over the unrenor- 
malized plane-wave impulse approximation, though the EMA predictions for differential 
cross sections still differ significantly from those of the full DWBA in heavy nuclei, even 
at high energies ||] . 

This paper discusses the application of these very convenient approximations and mod- 
ifications thereof to charged-current neutrino-nucleus scattering, for which the outgoing 
lepton is either an electron or muon. Because the neutrino is not charged and the weak 
interaction is short-ranged, the expansion in partial waves is less complicated and time- 
consuming than in electron scattering. But thus far the most significant observable in all 
the neutrino experiments cited above is the total cross section, and one might question 
whether the full DWBA is really required to obtain it. The two approximations just dis- 
cussed preserve the intuitive plane-wave picture of the scattering and are so much easier 
to apply that one really ought to retain them if possible. On the other hand, it is not clear 
that they will always work well; the EMA, for example, has been derived only for extremely 



relativistic particles [||, 10, 11 1, while the most energetic outgoing muons considered here 



are far from the relativistic limit. I show below however that although the Fermi function 
is accurate only for low-energy electrons in most nuclei, a simple modification allows the 
EMA to work remarkably well (better than in electron scattering), even in heavy nuclei 
and, for muons, even at very low energies. 



2 General Considerations 



The differential cross section for a neutrino of momentum k v to produce a muon or 
electron at scattered angle Q, and energy E, with the nucleus going from state \i) to state 
|/), can be written as a phase space factor times a sum of squared matrix elements, viz: 
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where G is the Fermi coupling constant, J is the charged weak nuclear current, vt is the 
4-spinor containing the momentum-space components of the left-handed neutrino plane 
wave with momentum k u , and if)Z (r) is the charged lepton (electron or muon) wave 

function with spin projection s and energy E = \fk 2 + M 2 in the electrostatic potential 

V = (3 - r 2 /R 2 ) r <R 

V = -Za/r r>R. (3) 

(R is the nuclear radius.) The wave function ip~ is an eigenstate of the full Hamiltonian, 
electrostatic potential included, that looks asymptotically like a plane plus an ingoing 
spherical wave. Expanded in partial waves, it has the form[|12|, [13] 
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where the quantum number j and the label I determine?] the the eigenvalue k = (/ — 
j)(2j + 1) of the operator K = —f3(2L -5 + 1) (J3 is the usual Dirac-equation matrix); 
^ijm(^) = Kx]&» i s a two-spinor with non-relativistic quantum numbers i, s = and j; 
and is the "inner phase shift" |0|. The radial wave functions P^i(r) and Q^(r) obey 
the equations 



~pZ(r)-[E + M-V(r)]Qf ;l (r) 
-Qf,(r)-[E-M-V(r)}pZ(r) . (5) 



Since the goal here is to test approximations that preserve the easily interpreted plane- 
wave formalism, I will use only the vector charge and not the full current in Eq. (|2|); the 
charge, for which full distorted-wave cross-section formulae can be displayed and calculated 
simply, should be sufficient. If the nuclear target states have SU(4) symmetry the major 
part of the axial-vector current will affect the matrix element in essentially the same way 



lr The total angular momentum j is a good quantum number but the orbital angular momentum / is not 
and serves as a label here. 



as the charge. There is no obvious reason to think that the rest of the current will change 
the results in a qualitative way. 

Using the upper and lower radial wave functions P^{r) and Qfi(r) from Eq. (||), one 
can define 



where 
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is the isovector transition transition density and the sum is over nucleons n. The angle- 
integrated cross section, which is of particular interest, takes a relatively simple form if the 
recoil of the nucleus is neglected, since k and E are then held fixed while f2 is integrated 
over, with the result that the cross section does not depend directly on the phase shifts 
and the partial waves do not interfere with one another. The general expression for a final 
state with arbitrary J 71 " (and initial state with + ), restricting the current to the vector 
charge, is 
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This expression will be used below for comparing exact results with approximations. I 
will take the density in Eq. (0) to have the form pji oc 5(r — R)YjM{r) for transitions 
from a + ground state to excited states with angular momentum J and projection M, 
so that, e.g., 1~ states correspond to the isobar analogs of Goldhaber- Teller collective- 
dipole excitations. There are several reasons for this choice: densities for most low-lying 
collective states are surface-peaked, placing the strength at the nuclear radius provides a 
severe test for the approximations below (which work best at small r) and, finally, densities 
for higher-energy noncollective states can be represented as superposed densities peaked at 
different points in the nucleus, so that a <5-function density should be sufficiently general. 



3 The Fermi Function 

It is convenient to think about the effects of an electrostatic potential on an outgoing 
particle in terms of a local effective energy and momentum inside the nucleus: 



E eS = E - V(0) , k eS = JE 2 ea - M2 (9) 



When k e ffR « 1, a Fermi function of Z and E that multiplies the outgoing wave is appro- 
priate. In this low "effective-momentum" limit, only s-waves contribute to the scattering, 
and they vary nearly linearly inside the nucleus; the Fermi function can then be taken to 
be the ratio of, e.g., the Coulomb wave to the free wave at the nuclear surface, a quantity 
given by 

F(Z,E) = 2(1 + ^)( 2 fcfl)-^) e ^ J^±^g , 70 = V /T3^2 , u = ±^. 

(10) 

The Fermi function is useful for low-energy electrons and is often employed in calcu- 
lations of beta decay, but is not likely to work well for muons except in light nuclei. The 
reason is that for nonrelativistic muons, k e R ~ y/2M^E c ff > y/2M^Vo, which even in 12 N, 
the product of a charge-exchange reaction on 12 C, implies that k e sR ~ .5 (large enough to 
cause a 10% error with the Fermi function above). Thus, even when the energy at infinity 
is so low that only s-wave muons are produced the usual Fermi function will generally not 
be not very accurate. Figure |l] shows the total cross section for exciting a (fictitious) + 
state in 208 Bi at 15 MeV with a surface-peaked transition density, as a function of outgoing 
lepton kinetic energy for both electrons and muons; the + multipole is the one with the 
largest s-wave contribution at low energies so the choice should maximize the accuracy of 
the Fermi function. Indeed, for electrons, the function does well for the lowest 5-10 MeV 
of electron energy. For muons, however, the approximation is never valid. Its accuracy 
can be improved significantly by using the s-wave solution from the potential in Eq. (|3|), 
rather than the pure Coulomb potential (see Ref. ||), but except at low energies where 
such functions are tabulated, that means solving the s-wave Dirac equation, a task that is 
simpler than carrying out the full DWBA but still relatively involved. 

One physical situation in which these issues are important is the scattering of elec- 
tron neutrinos produced by muon-decay. A few years ago an attempt was made to cali- 
brate an 127 I solar-neutrino detector by exposing it to a flux of muon-decay neutrinos at 



LAMPF|14j. The neutrinos have much higher energy than their solar counterparts and 
excite the nucleus in different ways, so that a good understanding of low- lying states in 
127 Xe is required to extract the solar-neutrino cross section. In addition, the outgoing 
electron often has too high an energy to allow use of the usual Fermi function (but not 
high enough for the EMA); if the transition density is surface-peaked and concentrated at 
an average of 6 MeV, the Fermi function gives values about 50% larger than the DWBA 
+ multipole cross section when only the vector charge is included. The difficulties with 
nuclear structure and electron distortion were discussed in ref. [15]; to treat the distortion 
for the + and 1 + components of the cross section, the authors relied on a fit [16 to ta- 
bles in Ref. Q, which use the potential in Eq. (J3j) rather than the pure Coulomb force to 
calculate the Fermi function. Here that procedure still misses the cross sections, though 
only by about 10 or 15%. The difference is due mainly to higher partial waves, which 
are affected differently by the electrostatic force than the s-wave that underlies the Fermi 
function. It would take a significantly larger difference, however, to alter the conclusion 
of Ref. [15 1 that the LAMPF cross sections are too large to be understood or interpreted. 



Figure 1: The total cross section for scattering from the ground states of 12 C and 208 Pb 
to (fictitious) + states at 15 MeV in 12 N and 208 Bi, with transition densities proportional 
to 5{r — R), as a function of outgoing lepton kinetic energy for both electrons (a) and 
muons (b). The solid lines represent the full DWBA results, the dotted lines neglect the 
Coulomb force completely, and the dashed lines are the Fermi-function approximations. 
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4 The EMA and a Modified EMA 

At higher energies, the Fermi function doesn't work even for electrons and we need 
a different approximation. In electron scattering, the EMA is used to untangle the elec- 
trostatic attraction of the electron to the nucleus from the single hard-photon exchange. 
As already mentioned, in its simplest form this approximation consists of shortening the 
electron wavelength k to k e g inside the nucleus, resulting in a larger effective momentum 
transfer, and rescaling the amplitude of the wave function by k e fi/k. The procedure's 
accuracy depends on scattering angle and decreases significantly as Z grows; it is not 
regarded as good in lead. But the application of the EMA to electrons produced by 
neutrino scattering, which is straightforward, should work better than its application in 
electron scattering. One reason is that only the outgoing wave experiences a Coulomb 
force; a more significant one is that the range of the weak interaction is short. In electron 
scattering regions far from the nucleus, where the local momentum is quite different from 
k e ff/k, affect the scattering matrix element because of long-range photon exchange. Here 
by contrast the only effects are from inside the nucleus. Figure ^ shows the total cross 
section for scattering from the ground states of 12 C and 208 Pb to Goldhaber- Teller (ana- 
log) resonances in 12 N and 208 Bi, along with the EMA predictions. The approximation is 
quite good above about 30 MeV in lead, leaving only a small kinematic window in which 
neither the EMA nor a Fermi function works. For electrons, therefore, one can accurately 
account for the effects of the electrostatic potential without using the DWBA except in 



restricted kinematic regions in heavy nuclei. 



Figure 2: The total cross section for scattering from the ground states of 12 C (a) and 
208 Pb (b) to (fictitious) 0+ states at 15 MeV in 12 N and 208 Bi, clS £1 function of outgoing 
kinetic energy and with transition densities proportional to 5(r — R). The solid lines are 
the full DWBA results, the dotted lines neglect the Coulomb force completely, and the 
dashed lines are the EMA results. 
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Muons are another story however; the best way to extend the EMA to massive particles 
is not immediately clear. The approximation and corrections to it have been derived from 
the high-energy limit of the Dirac equation, either through an eikonal approximation pj 
or an expansion in inverse powers of fcflnj. Both approaches begin by neglecting the mass 
of the electron and so separating the usually coupled left and right handed spinors, hence 
simplifying the equation. Here we are concerned with muon energies below a GeV, for 
which the rest mass is not negligible. The basic ideas behind the simplest approximation — 
that the Coulomb potential shortens the wavelength in the interaction region and increases 
the wave amplitude by focusing or defocusing the particles as they approach or leave — do 
not appear to hinge on the particles' rest mass, as long as the wavelength is short enough 
for a local momentum to have some meaning. Some version of the EMA should therefore 
apply to intermediate-energy muons as well as electrons. 

The shortening of the wavelength will clearly work the same way for the two kinds of 
particle, the only difference being the kinematic relation between k e g and E e g, but what 
about the amplitude rescaling? Should it still be k e g/k or some other factor that involves 
the muon mass? We can decide without extending the the full analysis of Refs. |], 10] 



by viewing the change in amplitude as follows: The use of a plane- wave approximation 
in the interaction region is equivalent to the assumption that the Coulomb potential does 
not cause the particles to alter their direction very much when they approach or leave 
the nucleus. It should therefore not strongly alter an outgoing muon wave packet, which 



asymptotically is spherical, after it leaves the nucleus except by slowing it down and 
thereby changing the average radial wavelength and amplitude as the wave moves to larger 
r. (In particular, it should not cause much reflection.) But the differential cross section 
is directly related to the radial integral of the square of the wave function over all r at a 
fixed solid angle (in the large time limit), and under the assumptions above the Coulomb 
interaction should not affect this quantity much once the outgoing wave moves beyond 
the nucleus. Since most of the effect it does have will be to redistribute flux from one 
angle to another, the total cross section should be nearly the same as if the wave kept the 
same radial wavelength it had inside the nucleus, i.e. as if the potential were equal to V(0) 
everywhere in space. For that situation the plane-wave approximation is exact to lowest 
order in Za if the effective momentum and energy k e g and E e g replace the real quantities 
k and E everywhere in Eq. @. In other words, besides shortening the wavelength in the 
matrix element A4, one should replace the phase space factor kE by k e gE e g Jl^]. In the 
Born approximation this is equivalent to a change in the muon wave function 



i.e. a rescaling of the amplitude by \Jk c gE c g/kE rather than k e g/k. I will call this ap- 
proximation the "modified EMA". 

With this alteration the EMA in fact generally works better for muons than for elec- 
trons at low energies for the same reason the Fermi function doesn't work as well: the 



muon mass guarantees that k e gR never drops below ■^2M^V{Qi)R ~ 0.5 in nitrogen and 
2.6 in bismuth. Figure ^ shows total cross sections for exciting several states with different 
angular momenta in Bi (again assuming a surface peaked transition density) alongside the 
modified EMA predictions and those of the usual EMA, in which the plane wave is scaled 
by keg/k. The usual EMA is pathological at low energies, but the modified EMA agrees 
extremely well with the exact solutions except for the + multipole, where it is off by a 
little over 5% at muon energies of 150 MeV. (There are fewer ways of coupling neutrino 
and muon angular momenta to + , which may accentuate errors in any one partial wave 
for that multipole.) The differences in lighter nuclei are less; in 56 Fe the agreement is 
almost exact for all multipoles except at very low energies, where the + multipole still is 
somewhat worse than the others. One might see similar behavior in the 1 + multipole, also 
"allowed", were the entire current included. But the modified EMA is still dramatically 
better than the usual plane-wave treatment, even for + transitions in heavy nuclei. 

Thus far the only significant laboratory experiments with muon neutrinos in this energy 
range are the LAMPF measurements of the cross section for neutrinos produced by the 
decay of pions in flight scattering from 12 C. Because Z is small, electrostatic effects are 
likewise small. Most of the difference between the DWBA and plane-wave approximations 



in this light nucleus was correctly accounted for in Ref. [18|, which rescaled the muon phase 
space as prescribed by the modified EMA without altering the nuclear matrix element 
(though the definition of k e g is reported incorrectly in the manuscript). Another situation 
worth examining, however, is the scattering of atmospheric muon and electron neutrinos 
from the iron in the SOUDAN proton-decay detector. Results of experiments there and 



Figure 3: Total cross sections for several multipole excitations (labeling the panels) from 
the ground state of 208 Pb to fictitious states at 15 MeV in 208 Bi, as a function of outgoing 
kinetic energy and with transition densities proportional to S(r — R)Yjm- The solid lines 
are the results of full DWBA, the dotted lines those with the Coulomb force neglected, 
the dashed lines those of the usual EMA used in electron scattering, and the dot-dashed 
lines those of the modified EMA. 
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at other places imply that the ratio of muon-neutrino flux to electron-neutrino flux in the 
atmosphere is about 50-80% of the expected ratio || . The analyses use the Fermi- gas model 
of the nucleus and don't generally account for final-state effects in quasielastic scattering. 
Table 1 shows the changes due to the electrostatic force (calculated in the modified EMA) 
from the simple Fermi-gas model predictions for the quasielastic production of muons, 
antimuons, electrons, and positrons in iron. I have chosen relatively low outgoing momenta 
and taken neutrino fluxes from Ref. [19|. The individual event rates actually change 
significantly for these momenta, but as is usually the case, the ratio of total /x to total e 
events hardly budges, so that one one apparently can safely neglect Coulomb effects even 
in a material as heavy as iron. 



Table 1: Calculated event rates for atmospheric neutrinos in 5 Fe, in arbitrary units, with 
and without Coulomb distortion of the final-state leptons. 
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5 Conclusions 

The chief result of this paper is that the effective-momentum approximation works 
better in neutrino scattering than in electron scattering, particularly for muons, where the 
approach has to be modified slightly to take into account the charged particle's mass. It 
probably doesn't work quite as well for differential cross sections but preliminary analysis 
(on + transitions with the vector charge) indicate that the modified EMA is still surpris- 
ingly good. The real scattering redistributes some flux from forward angles to places where 
the cross section is lower, but not as much as one might imagine. In any event, we are not 
likely to pay attention to details of weak-interaction differential cross sections in the near 
future. For total cross sections the modified EMA works admirably for muons down to 
low energies and the usual EMA just as well for electrons except in a region where a Fermi 
function is often adequate. In heavy nuclei there is a kinematic window for electrons in 
which neither approximation is completely sufficient, but it is small. Important examples 
from this window — neutrinos from muon decay on iodine is one — are not easy to find 
at present. 

This work was supported by the U.S. Department of Energy under grant DE-FG02- 
97ER41019. I thank E. Kolbe and P. Vogel for useful discussions and correspondence. 
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